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DYNAMIC  SNAP- THROUGH  OF  IMPERFECT 
VISCOELASTIC  SHALLOW  ARCHES1 


N.  C.  Huang2  and  W.  Nachbar3 


ABSTRACT 

Dynamic  snap-through  or  dynamic  buckling  of  imperfect 
viscoelastic  shallow  arches  with  hinged  ends  is  considered  under 
step  loads  of  infinite  duration.  Attention  is  principally  devoted  to  the 
influence  both  of  small  imperfections  and  of  small  amounts  of  damping, 
acting  together,  on  the  critical  loads.  For  the  problem  considered, 
the  Voigt  model  is  used  for  viscoelasticity,  the  deflection  is  represent¬ 
ed  by  the  first  two  harmonic  modes,  and  imperfections  have  the  shape 
of  the  second  (antisymmetric)  mode.  Results  obtained  by  numerical 
integration  of  the  differential  equations  show  that  the  critical  load 
exhibits  a  jump  discontinuity  in  the  limit  both  for  vanishing  imperfect¬ 
ion  and  for  vanishing  viscosity.  Critical  loads  for  slight  imperfect 
and  elastic  (inviscid)  arches  are  slightly  higher  than  those  from  the 
saddle  point  formula  of  Hoff  and  Bruce  (J.  Math.  Physics,  32.,  1954, 
276),  confirming  that  the  formula  gives  a  lower  bound  on  the  critical 
load.  However,  critical  loads  for  arches  with  slight  imperfection  and 
slight  viscosity  are  considerably  higher  than  for  the  elastic  arches. 
Another  closed-form  expression  is  shown  to  be  in  good  agreement  with 
these  results.  For  finite  amounts  of  viscosity,  the  critical  loads  tend 
rapidly  to  the  values  obtained  for  infinite  viscosity,  which  are  the  same 
as  the  critical  loads  for  quasi-static  buckling.  Apart  from  the  jump 
discontinuity  at  zero,  the  critical  load  for  any  viscosity  decreases  con¬ 
tinuously  and  monotonically  with  imperfection. 


1  Research  sponsored  by  the  Air  Force  Office  of  Scientific  Research, 
Office  of  Aerospace  Research,  United  States  Air  Force  under 
AFOSR  Grant  AF-AFOSR  1226-67. 
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California,  San  Diego,  La  Jolla,  California 
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perfect  arch  dynamic  snap>through  load,  Eq.  (32) 
quasi-static  snap-through  -  see  last  column  of  Table  3 
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potential  energy,  E4.  (24) 
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2  bx  (t)  ;  first  harmonic  of  deflection 
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=  z  (0) :  imperfection  amplitude 
value  of  z  at  critical  point 

approximate  shortening  of  arch  length,  Eq.  (5) 
effective  viscosity  for  Voigt  solid 
nondimens ional  viscosity,  Eq.  (15) 
mass  of  arch  per  unit  length 
nondime  ns  ional  time,  Eq.  (15) 
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1.  INTRODUCTION 


Iri  this  paper,  the  dynamic  snap-through  or  dynamic 
buckling  behavior  of  imperfect  viscoelastic  shallow  arches  under  a 
suddenly  applied  load  of  infinite  duration  is  considered.  The  critical 
load  is  defined  as  the  smallest  magnitude  of  load  under  which  the  arch 
can  deform  to  a  snapped-through  position. 

The  problem  of  the  dynamic  buckling  of  sinusoidal  elastic 
shallow  arches  has  been  studied  by  Hoff  and  Bruce  [l]  4  and  by 
Simitses  [2J  with  analyses  based  upon  the  behavior  of  the  trajectories 
describing  the  motion  on  the  potential  surface.  Since  small  geometric 
imperfections,  or  small  disturbing  loads  or  initial  velocities,  affect 
the  shape  of  the  potential  surface  only  to  the  second  order  in  small 
quantities,  it  is  possible  to  include  the  effects  of  such  perturbations  by 
an  analysis  based  upon  the  potential  surface  for  the  perfect  structure. 
This  analysis  for  the  arch  shows  that,  if  the  arch  is  not  very  flat,  the 
critical  load  of  a  perfect  arch  having  the  shape  of  sin  x  will  be  assoc* 
iated  with  an  asymmetric  buckling  mode  shape  sin  2x.  This  critical 
load  is  obtained  from  the  condition  that  it  is  possible  for  i  trajectory 


Numbers  in  brackets  designate  references  at  the  end  of  the  paper. 


to  reach  a  saddle  point  on  the  potential  surface.  The  critical  load  is 
that  load  for  which  the  potential  energy  at  the  saddle  point  becomes 
equal  to  the  total  energy. 

Several  reasonable  objections  have  been  raised  to  this 
potential  surface  method  of  analysis.  One  objection  is  that  the  method 
is  limited  to  problems  with  a  small  number  of  degrees  of  freedom. 

The  arguments  used  in  the  potential  surface  method  are  essentially 
geometric,  and  it  is  much  more  difficult  to  deal  with  the  hyper¬ 
surfaces  that  are  encountered  if  more  than  two  degrees  of  freedom  are 
considered.  We  will  not  attempt  in  this  paper  to  generalize  this  method 
of  analysis  but  will,  in  fact,confine  our  discussion  to  a  specific  system 
with  two  degrees  of  freedom  where  the  method  has  already  proven 
feasible. 

Another  objection,  and  the  one  to  which  the  major  attention 
of  this  paper  will  be  devoted,  is  that  the  saddle  point  criteria  does  not 
assure  that  the  solution  to  the  differential  equations,  with  its  initial 
conditions,  does  indeed  have  a  trajectory  on  the  potential  energy  surface 
which  goes  into  the  saddle  point  at  the  critical  load.  Hoff  and  Bruce 
pointed  to  this  deficiency  in  the  following  statement  of  their  paper  re¬ 
garding  the  criticsl  value  of  their  load  parameter  Q,  ([l],  pp»  284-5): 
"If  Q  has  this  value  one  can  say  that  it  is  possible  for  the  initially 
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disturbed  system  to  be  displaced  to  the  stable  'buckled'  equilibrium 
configuration,  (0,  ezj,  if  it  follows  the  appropriate  path  on  the  energy 
surface.  "  We  have  underlined  the  words  possible  and  if.  The  signifi¬ 
cance  of  these  words  is  amplified  in  the  work  done  in  [2],  where  the 
Hoff  and  Bruce  critical  load  is  not  regarded  as  a  rigorous  critical  load 
but  as  a  lower  bound  or  minimum  possible  critical  load  for  dynamic 
snap-through  of  the  arch. 

The  question  then  arises:  if  the  saddle  point  critical  load 
is  a  lower  bound  on  the  true  critical  load,  how  good  an  approximation 
to  the  critical  load  is  it?  To  investigate  this  question,  it  is  necessary 
to  integrate  the  governing  differential  equations,  to  look  for  the  minim¬ 
um  load  at  which  the  solution  exhibits  snap-through,  and  to  compare 
this  minimum  load  with  the  predictions  of  the  potential  surface  analysis. 
In  carrying  out  this  program,  we  cannot  confine  ourselves  to  perfect 
arches  of  sinusoidal  shape.  Indeed,  it  is  readily  seen  that  if  the 
differential  equations  and  subsidiary  conditions  are  symmetric,  then 
the  solutions  obtained  will  be  symmetric,  and  the  trajectories  of  these 
solutions  could  never  approach  a  saddle  point  that  is  associated  with  an 
unsymmetrical  mode.  It  is  necessary  then  to  consider  imperfections 
or,  at  the  least,  to  consider  the  limiting  behavior  of  an  imperfect 
arch  as  the  amplitude  of  the  imperfections  tends  towards  zero.  A  more 
elaborate  analysis  [2]  has  already  indicated  that  second  harmonic  (or 


first  antisymmetric)  imperfections  will  have  the  greatest  effect  on  the 
critical  load  and  so,  in  the  present  paper,  only  the  first  and  second 
harmonics  are  considered  to  describe  the  initial  shape  and  subsequent: 
deformation  of  the  arch. 

Viscoelasticity  is  introduced  here  because  it  is  a  more 
realistic  material  description  for  dynamic  problems  than  the  assumption 
of  ideal  elasticity,  and  also  because  it  often  makes  the  analysis  less 
ambiguous,  as  is  illustrated  by  our  study  of  dynamic  buckling  of  a 
single  degree-of-freedom  structure  [3]  .  Perhaps  more  important, 
however,  is  that  we  want  to  see  whether  the  critical  loads  obtained  for 
the  viscoelastic  material  in  the  limit  as  the  viscosity  parameter  tends 
to  zero  are  the  same  as  the  critical  loads  obtained  under  the  same  con¬ 
ditions  for  the  elastic  material  in  which  the  viscosity  is  equal  to  zero. 

In  the  case  of  imperfections  we  already  know  that  the  limit  of  the 
critical  load  for  vanishing  imperfections  may  not  be  equal  to  the 
critical  load  for  zero  imperfections.  We  then  may  have  a  jump  between 
the  limiting  value  and  the  value  at  zero.  We  want  to  see  whether  such 
a  jump  is  exhibited  also  with  respect  to  the  viscosity  parameter. 

We  will  also  present,  in  addition  to  the  investigations 
already  described,  results  on  the  effects  of  finite  imperfections  and 
finite  viscosity  on  the  critical  load.  In  this  paper,  we  will  discuss  all 
these  questions  only  with  regard  to  the  following  problem :  the 


arch  is  considered  hinged  at  both  ends;  the  loading  is  sinusoidal  in 
shape,  and  is  of  constant  magnitude  and  infinite  duration  in  time;  the 
material  of  the  arch  is  considered  to  be  a  Voigt  solid. 

2.  BASIC  EQUATIONS 

Consider  a  shallow  arch  of  span  L,  as  shown  in  Fig.  1. 

The  cross  section  of  the  arch  is  assumed  to  be  uniform  along  the  span, 

with  cross-sectional  area  A  and  moment  of  inertia  I  about  the  central 

axis.  The  load  q(x)  per  unit  length  is  suddenly  applied  at  zero  time 

and  maintained  at  constant  magnitude  thereafter.  At  time  t,  the  central 

axis  of  the  arch,  having  original  height  wq(x),  is  deformed  to  a  height 

w  (x,  t) .  With  use  of  d'Alembert's  principle,  the  dynamic  load  is 
d3  w 

q(x)  +  p  ^  '"z"  positive  in  the  -  w  direction,  as  shown  in  Fig.  2,  where 
p  is  the  mass  of  the  arch  per  unit  length.  If  we  denote  the  horizontal 
thrust,  vertical  shear  and  bending  moment  by  H,  V  and  M,  respect¬ 
ively,  the  equations  of  equilibrium  can  be  obtained  directly  from  Fig.  2: 

V'  =  -  (q  +  pw)  (1) 

and 

M'  =  V  -  Hw'  ,  (2) 

where  (  )'  =  (  )  and  f)  =  ^  (  )  • 

If  the  arch  is  elastic,  the  Euler-Bernoulli  assumption  gives 

M  =  EI(w"  -  w  ")  (3) 

o 
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where  E  is  Young's  modulus.  From  Eqs.  (1),  (2)  and  (3)  the  follow¬ 
ing  equation  is  obtained: 

E I  ( v /"  -  v/J"'  )  +  H  w"  +  q  +  p  w  =  0  .  (4) 

For  shallow  arches,  the  shortening  of  the  length  of  the 
arch  during  deformation  can  be  written  approximately  as 

A  -  2  J  [(*.')  -  ("')  ]d*  ■  |5> 

o 
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w{x,  t) 


(11) 


I  v1/3  ® 

-2\t)  s,  b»<*>8i" 

n=l 


ntr  x 
L 


where  the  deflections  have  been  normalized  with  respect  to  the 

radius  of  gyration  (I/A)1/**,  so  that  the  constants  a  and  the  functions 

n 

b  (t)  are  dimensionless.  Further,  let  the  load  intensity  be 
n 


represented  by 


q(x)  =  E  f  sin  .  (12) 

.  n  L 

n=l 


Substituting  Eqs.  (10),  (11)  and  (12)  into  Eq.  (7),  we  obtain 

n*  E(bn  -  aj  +  n*  b^E  ^  If  <b.*  -  >.“)]  ♦  [(4)*\  +  2Pb,]=  0 


for  n  =  1,  2,  ... 


(13) 


We  shall  use  the  Voigt  model  to  represent  the  viscoelastic 

property  of  the  arch.  The  elastic  constant  E  in  Eq.  (13)  is  therefore 

d 

replaced  uy  the  differential  operator  E^  +  V  >  where  E^  is  the 
effective  spring  constant  and  TJ  is  the  effective  viscosity  of  the 
material.  Consequently,  Eq.  (13)  becomes 

"‘(Eo+”lr)  K  -  \)  *  a‘\  (Eo  +  (bJ*  -  aj3)] 

+  [(t)  fn  +  2-b„]  -  0 

for  n  -  1,  2,  .  .  .  .  (14) 

Let  us  introduce  the  following  dimensionless  quantities 
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and  notation: 


w 


■&■ 


=  'S1\ 

*->_  -si 


R  -  L‘  rA\1/a£ 

n  ‘  2tt*E  I  \  I  )  n  ’ 

o 

a 

ir 

T  =  T? 

(tf‘  •' 

IT*  (  I  \Va 

*  '  ”  s?  (r?)  and 

C)  ■ 

£<>  • 

(15) 

Equation  (14)  becomes 
d 


n‘  (* +  *  w)  K  -  *.» +  *' \  ( 1  +  *  h)  *  Ds  (V  - •,*)] + 


+  R  +  b  =0 
n  n 


for  n  =  1,  2 .  (16) 

We  shall  assume  that  the  load  distribution  is  sinusoidal 


Solutions  to  these  equations  will  be  sought  for  the  initial  conditions: 

y(°)  =  yo>  z{0)  =  *o,  y  (0)  =  yo  and  z(0)  =  zq  .  (19) 

We  can,  without  restricting  the  essential  problem,  consider  y  >  0, 

o 

z  5  0  and  r  5  0 . 
o 

The  critical  points  of  Eqs.  (17)  and  (18)  are  obtained  by- 
setting  y^y=z=z=0  in  these  equations.  The  critical  points  defined  by 
the  pairs  of  values  (y^,  z^)  satisfy  the  following  equations  simultan¬ 


eous  ly: 


y  -  y  +  yfy2-y2  +  4z2-4za^  +  r  =  0  , 

e  o  e  \  e  o  e  o  / 

4fz-z)+zfy3-y2  +  4z3-4za^=0 
\  e  o  /  eV'e  o  e  o  / 


(20) 


(21) 


After  z^  is  eliminated  from  Eqs.  (20)  and  (21),  we  obtain 

(3  ye  +  yo-  r)  [r  +  »e  *  yo  +  r«  (  y.“  -  yo’  *  4,o*)]  +  “*0“  y.*  =  °' 

(22) 

The  energy  equation  can  be  derived  from  Eqs.  (17), 

(18)  and  (19)  by  integration  and  addition.  It  can  be  written  as 

U  +  K  +  D  =  0  .  (23) 

where 

U  =  potential  energy 

=  \  (y  -  y Qf  +  8(z  -  zQ)a  +  ^  JV  -  Yq  +  4(za  -  zQa)J  + 


+  r  (y  -  yo)  ; 


(24) 
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$¥**»>»»> .  ...  .  ,  .  <i  .  ,■  ' *• 


K  =  kinetic  energy 


*-  '  -Sv 


-£(*■  -  V)  +  I(i3-  *.*)  * 


D  =  total  energy  dissipation 


L 

=  X  J  [y2  +  l6za  +  2(yy  +  4zz)a]  dt 


Note  that  the  datum  of  the  potential  surface  U  has  been  chosen  so 
that  the  initial  value  of  U  is  zero.  Equations  (20)  and  (21)  defining 

the  critical  points  of  Eqs.  (17)  and  (18)  are  the  same  as  the  equations 

dU  dU 

=  0  and  =  0 ,  respectively. 

3.  PERFECT  ARCHES  WITH  PERFECT  INITIAL  CONDITIONS 

(*  =  *_  =  0) 
o  o _ 

From  Eqs.  (17),  (18)  and  (19),  we  find  that  unsymmetric 

deformation  of  the  arch  can  be  caused  by  small  antisymmetric  initial 

velocities  or  imperfections  in  geometry.  On  the  other  hand,  if  the 

initial  conditions  are  entirely  symmetric,  i.  e. ,  if  z  =  z  =0,  we 

o  o 

can  conclude  from  Eq.  (18)  that  z  (t)  =  0.  Therefore,  under  the 
constraints  assumed  above,  the  deformation  of  a  perfect  sinusoidal 
arch  under  perfect  initial  conditions  remains  symmetric,  and  the  motion 
of  the  structure  has  only  one  degree  of  freedom  which  is  described  by 
y(t).  The  differential  equation  and  initial  conditions  are  taken  as 


1 


W-  yM t 
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(27) 


y  +  x (l  +  2ya )  y  +  y  -  yQ  +  y3  -  yoJy  +  r  =  0  .  (27) 

y  (0)  =  yQ  and  y(0)  =  0  .  (28) 

The  critical  points  in  this  problem  are  y  =  0,  y  =  y  , 
where  y^  is  a  root  of  the  following  cubic  equation  (corresponding  to 
the  roots  of  the  square  bracketed  term  in  Eq.  (22)  when  =  0) : 

ye3  +  (1  “  O  ye  +  r  ‘  yo  =  0  *  (29) 

Typical  solutions  of  Eqs.  (27)  and  (28)  are  shown  in 
Fig.  3.  These  solutions  were  obtained  by  numerical  integration  (see 
part  (b)  below).  For  positive  X,  these  curves  oscillate  about  one  of 
the  stable  equilibrium  points  (which  are  spiral  points  in  the  phase  plane) 
with  decreasing  amplitudes  because  of  continuous  energy  dissipation.6 
When  the  load  parameter  r  passes  through  a  critical  value8  r^  there 
occurs  a  sudden  diversion  of  these  curves  from  the  neighborhood  of 
one  equilibrium  point  to  the  neighborhood  of  the  other. 

For  problems  described  by  a  single  variable,  the  analysis 
of  dynamic  snap-through  is  facilitated  by  the  use  of  potential  curves. 


6  When  the  value  of  x  is  very  small,  the  structure  may  snap  back 
after  the  snap-through  (see  [3]  ). 

8  In  order  to  be  able  to  distinguish  the  critical  load  values  obtained 
by  numerical  integration  of  the  differential  equations  as  opposed  to 
closed-form  expressions  for  r  obtained  for  particular  cases  by 
use  of  the  potential  surface,  we  will  adopt  the  convention  of 
identifying  the  latter  by  special  subscripts  (e.  g.  rcn,  r  ,  r  , 
r  ),  and  we  will  reserve  the  notation  r  for  the  former. 


i 


I 


- 


Equation  (24)  becomes 

U  =  J  (y  -  yQ)a  +  ^  (ya  -  yQa)a  +  r(y  -  yQ)  •  (30) 

Using  r  as  a  parameter,  we  can  plot  U  as  a  function  of  y,  as 
shown  in  Fig.  4.  A  discussion  of  the  application  of  the  potential 
curver  for  dynamic  buckling  analysis  in  a  dissipative  system  with  a 
single  degree  of  freedom  is  given  in  [3]  .  In  the  following  parallel 
discussion,  snap-through,  as  a  function  of  X ,  is  discussed  as  three 
cases. 


w. 


(a)  Elastic  case  (X  =  0). 

No  solution  of  Eqs.  (27)  and  (28)  can  produce  a  positive  U, 

in  view  of  Eq.  (23).  With  r  increasing  from  zero,  then,  the 

minimum  of  the  potential  curve  on  the  left  of  the  U  axis  (Fig.  4) 

will  not  be  reached  by  a  solution  until  the  local  maximum  value  of 

U  to  the  right  of  this  minimum  satisfies  the  condition  U  £  0 . 

e 

Therefore,  the  critical  load  can  be  obtained  by  setting  the  value  of 
U  equal  to  zero  at  the  equilibrium  point  at  which  this  maximum 


is  attained.  This  gives,  from  Eq.  (30), 

ye)[2  +  X  yo)  ] 


r  = 


By  eliminating  y^  between  Eqs.  (29)  and  (31),  we  obtain  an 


expression  for  the  critical  load  r 


cp 


=  W  [2yo  (9  +  X  *  (2yo“  *  3)  (4yo“  -  6)V  ]  • 


cp 
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(31) 


se 

i 


The  value*  of  r  are  plotted  against  yQ  by  one  of  the  solid 
lines  in  Fig.  5. 

From  Eqs.  (29)  and  (30),  the  potential  energy  U  at  all 
critical  points  is  found  to  be 

ue  ■  -  (t.  -  O'  [3  (t.  +  T  y0  )s  + !  (3  •  2y0* )]  •  (33) 

If  y  is  less  than  ~  (6)^3  =  1.224,  Eq.  (33)  shows  that  U  is 

O  4  6 

negative  at  all  critical  points  ;  hence,  the  arch  always  deforms 
to  the  inside -out  position  without  snap-through. 

(b)  Viscoelastic  case  (0  <  x  <  ® ) . 

For  the  viscoelastic  case,  it  is  necessary  to  use  numerical 
methods  to  solve  the  differential  equation,  Eq.  (27),  with  the  initial 
condition,  Eq.  (28).  A  Runge-Kutta-Gill  method  [4]  has  been  used 
for  this  purpose.  The  critical  load  rc  is  determined  by  the  con¬ 
dition  that  the  arch  begins  at  that  load  to  deform  to  the  snapped- 
through  position.  The  results  for  cases  X  =  0. 1  and  0.2  are 
shown  in  Fig.  5  by  solid  lines. 

A  certain  amount  of  energy  is  dissipated  before  the  struc¬ 
ture  reaches  the  point  of  snap-through.  Hence,  the  critical  load 
for  the  viscoelastic  case  is  higher  than  that  in  the  corresponding 
elastic  case.  In  the  viscoelastic  case,  there  is  also  a  lower  limit 
for  y^  such  that,  for  any  value  of  y^  less  than  this  lower  limit. 
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the  arch  deforms  to  the  inside-out  position  without  snap-through, 
(c)  Quasi-static  case  (X  =  a). 


'..’tv  ;-rr. 

'  vSr  j. 


The  quasi-static  case  was  first  studied  by  Fung  and 

Kaplan  [5].  The  same  results  are  obtained  here  in  the  limit  of 

infinite  viscosity.  When  X  =  ® ,  snap-through  can  occur  only 

when  the  maximum  of  the  potential  curve  (Eq.  30)  becomes  a 

daU 

point  of  inflection,  i.  e.  ,  at  y  =  y^,  =  0  .  From  this  con¬ 


dition  we  have 


The  quasi-static  critical  load  r  can  then  be  derived  from 

cs 


Eqs.  (29)  and  (34): 


yo  *  |  Pftyo*  -  1  ) 


la  Fig.  5,  r  is  plotted  against  y^  by  a  solid  line.  Note 
that  when  yQ  <  if  there  cannot  be  a  point  of  inflection  on  the 
potential  curve,  Hence,  for  y  <  1,  the  arch  will  deform  to 
an  inside-out  position  under  any  positive  load. 


4.  ARCHES  WITH  INFINITESIMAL  GEOMETRIC  IMPERFECTION 

(*  =  0  +  ) 

_ o _ 

The  perfect  arch  discussed  in  the  last  section  is  an  overly 
idealized  model.  It  is  known  that  the  critical  load  may  change 


drastically  due  to  the  effect  of  small  geometric  imperfections. 
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General  studies  of  this  behavior  for  elastic  structures  have  been  made 
by  Koiter  [6]  and  Budiansky  and  Hutchinson  [7  ].  For  our  study,  we 
shall  again  assume  that  y^  =  =  0  and  consider  the  following 

three  cases. 

(a)  Elastic  case  (x  =  0) . 

The  effect  of  infinitesimal  geometric  imperfections  for 
the  problem  of  elastic  shallow  arches  was  first  studied  by  Hoff 
and  Bruce  Cl 3  with  the  use  of  the  potential  surface  for  the  perfect 
arch.  Their  results  are  to  be  obtained  here  through  consideration 
of  the  limiting  behavior  of  an  imperfect  arch  with  vanishing 
imperfections. 

When  the  imperfection  is  small,  so  that  <<  y^a , 

the  critical  points  (y^,  z^)  can  be  determined  by  the  following 

two  simultaneous  equations  derived  from  Eqs.  (20)  and  (21) : 

y  -  y  +  y(ya-y3  +  4z3)+r=0  ,  (36a) 

e  o  e\e  o  e  / 

4  (ze  '  ■‘o)  +  %  (*.’  -  +  )  =  °  •  (36b> 

Consider  first  that  y  is  some  fixed  number  and  that 

o 

z  tends  to  zero.  In  the  limit  as  z  tends  to  zero,  z  either 
o  o  e 

approaches  zero  or  approaches  a  nonzero  limit  value.  If  z 

e 

approaches  zero,  then  Eq.  (36a)  reduces  to  Eq.  (29).  If  z 

e 

approaches  a  nonzero  limit,  however,  then  Eq.  (36b)  becomes 
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Ze  (».*  *  V  +  4*.‘  +  4  )  *  0  •  <37> 

When  s  2,  we  see  that  the  expression  in  parenthesis  in 

Eq.  (37)  is  positive.  Equation  (37)  would  then  imply  z  =  0. 

Hence,  the  limit  of  z  as  z  tends  to  zero  can  be  different 

e  o 

from  zero  only  if  y^  >  2. 

Equations  (36a,  b)  can  be  combined  to  give 

z  (3y  +  y  •  r  )  •  4y  *  =0 

e  \  7e  ’o  )  7e  o 

If  the  limit  of  z  as  z  tends  to  zero  is  not  zero,  this 
e  o 

equation  gives  in  the  limit 

y.  *  ?<r  -  V  •  y„  *  2  <38> 

and  it  then  follows  from  Eq.  (37)  that 

z.  -  *  I  [yo*  -  4  -  ?  <r  -  V*  T  ■  yo  -  2  •  (39» 

Of  the  three  critical  points  on  the  y-axis,  two  are  local 
minima  and  one  is  a  local  maxim  am.  We  shall  denote  the  local 
minima  by  Mx  and  Mg,  where  Mx  is  the  minimum  closest  to  the 
initial  point.  Equations  (38)  and  (39)  define  two  saddle  points. 


svmmetric  about  z  =  0. 


From  Eq.  (24),  the  potential  energy  for  small  is 


U  '  2  (y-yo)  +  ^  *i{y‘-yo  *  *Z‘)  *  r( 


Y  -  ro  )  •  (40) 


The  potential  energy  U  =  Ug  at  the  saddle  points  is  found  from 
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5* 

f 


Eqs.  (38),  (39)  and  (40) : 

ue  =  i[(r-  4yc)‘- 24]  •  y0£  2  •  <41> 

Under  small  load,  the  trajectory  of  the  motion  on 
(y,  z)  plane  would  move  in  the  neighborhood  of  .  The  neigh¬ 
borhood  of  M„  is  not  accessible  until  the  load  is  increased 
sufficiently  so  that  there  will  exist  a  trajectory  leading  to  the 
neighborhood  of  and  along  which  U  is  nonpositive.  A 
necessary  condition  for  the  existence  at  such  a  trajectory  is 
given  by  consideration  of  the  saddle  point  on  the  potential  stir- 
face.  The  considered  trajectory  is  seen  to  be  impossible 
unless  r  is  large  enough  so  that  the  saddle  point  value  of  U 

e 

of  Eq.  (41)  is  nonpositive.  Hence,  the  critical  condition  of 

Hoff  and  Bruce  for  snap-through  is 

Ue  =  0  ,  (42) 

at  the  saddle  points,  and  the  associated  critical  load  is  found 

to  be  r  : 
cm 

r  =  4y  -  2  (6)^2  ,  y  a  2.04  (43) 

cm  o  7  o 

When  Eq.  (43)  is  used  to  substitute  for  r  in  Eqs.  (38) 

and  (39),  the  coordinates  of  the  saddle  points  at  r  =  r  are 

cm 

found  to  be 

ye  =  yQ  -  |  (6)1/3  (44a) 


17 


(44b) 


Equation  (44b)  shows  that  the  saddle  point  does  not  exist  if 

is  less  than  (5/(6)^S)  =  2.04.  At  yQ  =  5/(6)^2,  Eqs.  (44a, b) 

give  z  =0,  y  =  l/^)1/3  at  r  =  8/(6)^^ .  This  point  is  also  a 
6  e  e  cm 

root  of  Eq.  (29).  It  is  the  critical  point  which  corresponds  to  the 

local  maximum  of  U,  as  is  discussed  in  case  (a)  of  Section  3. 

Moreover,  the  value  of  U  at  that  local  maximum  is  zero.  Hence, 

at  y  =  5/(6)1/3,  we  have  r  =  r  =  8/(6)3^2,  the  value  of  r 
7  cm  cp  cp 

being  computed  from  Eq.  (32). 

Consequently,  Eq.  (43)  holds  only  for  y^  >  2.04.  For 

y^  £  2.04,  the  critical  load  for  the  slightly  imperfect  arch  is  the 

same  as  that  for  the  perfect  arch;  it  is  r  as  given  by  Eq.  (32). 

In  Fig.  7,  the  result  of  Eq.  (43)  is  shown  by  one  of  the 

dotted  lines.  It  is  seen  that,  if  r  is  an  accurate  measure  of 

cm 

critical  load,  then  for  yQ  >  2.04  the  critical  loads  drop  consider¬ 
ably  due  to  the  effect  of  small  antis ymmetrical  imperfections. 

Since  the  above  method  of  determining  the  critical  load  of 
an  imperfect  arch  does  not  examine  the  actual  motion  of  the 
structure,  it  does  not  enable  one  to  assert  that  the  trajectory  will 
pass  through  the  saddle  point  even  if  Eq.  (42)  is  satisfied.  Equa¬ 
tion  (42)  is  a  necessary  condition  for  snap-through,  but  it  io  not  a 


sufficient  condition. 


In  order  to  clarify  this  point,  we  have  obtained  solutions 
to  Eqs.  (17),  (18)  and  (19)  numerically  by  using  very  small 
values  for  the  geometrical  imperfection  parameter  z^  in  the 

_  3  .  Q 

range  10  to  10  .  The  critical  loads  obtained  from  the 

.  4 

numerical  solution  for  z  =10  and  various  y  ,  as  well  as 
those  from  Eqs.  (32)  and  (43),  are  shown  in  Table  1  and  Fig.  6. 

It  is  found  that  the  calculated  r^  is  in  agr  eement  with  the  pre¬ 
dictions  of  Eq.  (43)  to  within  an  error  not  exceeding  a  few  percent. 
These  results  confirm  that  it  is  impossible  to  produce 

snap- through  if  the  load  is  lower  than  the  value  r  given  by 

cm 

Eq.  (43).  Consequently,  the  result  of  Hoff  and  Bruce  furnishes 
a  lower  bound  and  a  good  approximation  to  the  critical  load  for 
elastic  arches  with  small  imperfections,  as  was  implied  in  [2].  The 
following  observations  were  made  from  the  numerical  results: 

(i)  the  arch  may  oscillate  nonlinearly  for  a  period  of  time  before 
it  snaps  through;  (ii)  after  snap-through,  the  arch  may  snap  back 
to  its  original  shape. 

If,  for  a  fixed  y  >  2.  04,  the  calculated  critical  load  r 
o  c 

is  sought  for  very  small  values  of  z^,  anomalous  results  maybe 

obtained.  Typical  of  such  results  are  those  shown  in  Table  2  for 

y  =  4.  5.  Values  of  r  appear  to  increase  as  z  gets  smaller, 
o  c  o 

The  reason  for  this  is  that  the  potential  surface  over  which  the 
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* 


=*  * 


trajectory  travels  is  very  "flat".  This  follows  because,  for  small 

z  ,  at  z  =  z  is  0(z  ).  Hence,  in  the  absence  of  damping,  the 
o  oz  o  o 

arch  can  oscillate  for  very  long  periods  before  snap-through.  If 
the  calculations  on  the  computer  are  not  continued  long  enough 


for  these  very  small  values  of  z^,  then  we  run  the  risk  of  not 
discerning  motions  which  are  actually  buckled  motions,  simply 


because  they  are  not  observed  for  long  enough  time. 


(b)  Viscoelastic  case  {0  <  X  <  00 ) 

Values  of  obtained  by  numerical  solution  of  Eqs.  (17), 

(18)  and  (19)  with  X  =  0.  1  and  X  =  0.  2  are  shown  in  Fig.  7 

_  4 

by  dotted  lines  for  very  small  imperfections  (z  =10  was  the 

o 

value  used)  and  by  solid  lines  for  zero  imperfection.  The  closed- 
form  expressions  obtained  for  the  perfect  elastic  arch  (Eqs.  (32) 
and  (35))  are  likewise  represented  by  solid  lines,  while  those  for 
the  infinitesimally  imperfect  elastic  arch  (Eq.  (43)  and  Eq.  (48) 
below)  are  represented  by  dotted  lines.  As  was  noted  in  (c)  of 
Section  3,  the  limit  x  =  »  gives  the  same  results  as  for  quasi¬ 
static  loading  of  r  n  elastic  arch. 

These  curves  show  that  for  finite  values  of  X,  and  for 
values  of  y^  above  a  certain  lower  limit,  r^  drops  markedly  with 
the  introduction  of  a  small  imperfection.  The  lower  limit  for  y  , 
which  is  2.  04  for  the  elastic  arch,  is  seen  to  be  higher  for  nonzero  X. 


Moreover,  the  furious  behavior  exhibited  in  Table  2  for  the 

elastic  arch  with  very  small  values  of  is  not  shown  for  the 

viscoelastic  arch;  in  the  latter  case,  the  limit  value  of  r  is 

c 

approached  continuously  in  the  calculations  as  goes  to  zero. 

(cf.  also  Fig.  8). 

For  values  of  y  not  much  above  the  lower  limit,  the 
o 

dotted  curves  for  X  =  0.  1  and  x  =  0.2  in  Fig.  7  are  seen  to 

merge  rapidly  into  a  calculated  dotted  curve  for  X  =  ®  in  Fig.  7. 

Other  data  for  smaller  values  of  X,  not  shown  in  the  figure, 

exhibit  the  same  behavior.  The  data  for  the  elastic  arch,  however, 

follow  close  to  the  curve  r  for  all  y  >  2.  04.  The  transition 

cm  o 

between  these  two  types  of  behavior  occurs  at  values  of  X  ^  10  , 

but  computations  are  again  difficult  for  this  regime. 

(c)  Quasi-static  case  (X  =  ®) 

Fung  and  Kaplan  (  [5],  pp.  14  and  15)  considered  in  some 

detail  the  snap-through  of  shallow  arches  due  to  antisymmetric 

bifurcation  in  quasi-static  loading  if  y  is  greater  than 
’  1/2 

—  (22p  =  2.  34.  To  render  the  present  paper  reasonably  complete, 

we  shall  derive  their  results  from  our  equations  following  the 
approach  given  in  [2]  . 

For  quasi-static  equilibrium  at  a  point  (y^,  z^),  under  a 
load  r,  Eqs.  (20)  and  (21)  must  be  satisfied  and  in  addition  the 
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potential  energy  U  must  obey  the  following  stability  condition  there: 


d3U  daU 

ay3 


a~u  /  Sa  U  \ 

dz3  \  dy  dz  / 


>  0 


(45) 


For  zq  vanishingly  small,  we  find  from  Eq.  (24),  the  following  partial 
derivatives  at  points  with  z  =  0  • 

S3U  2  a 

U  =  T-2  =  1  +  3ya  -  y  3  . 

yy  oy*  7  7o 

A  r  3  2  ^ 

U  =  t~2  =  4Ly-ya  +  4]  . 

ZZ  OZz  o 

tt  -  a2u 

U  —  ■%  -  0 


zy  dydz 


The  condition  (45)  can  then  be  written  as 


.*■ 

i 


& 

^  |? 


(l  +  3ys  -  y‘  )  (y3  ”  Y q  +  4 )  >  0 


u 


yy 


u 


ZZ 


(46) 


If  r  is  sufficiently  small,  so  that  y  -  y  is  small, 

o  e 

the  inequality  (46)  is  evidently  satisfied.  As  r  increases  and  y 

e 

decreases,  either  or  Uzz  will  become  negative,  indicating  a 

loss  of  stability.  The  critical  values  of  y  are  called  y  and  y 

©  C  1  C  Cm  ! 


U  =  0  : 

yy 


u  =  0  : 

ZZ 


_  /yo 


"cl 


rc2 


-  1  1/3 

\  ) 


■  iVo  ■  4) 


1/9 


(47a) 


(47b) 


It  is  easily  verified  that  at  y  =  ^  (22)1/3,  we  have  y  ,  =  y  ,,  and 

o  2  7cl  7c2 


mi- 
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that  <  Yc|  for  VQ  <  (22)1/3 .  Consequently,  the  critical 

value  of  y  is  y  ,  if  y  <  ^r(ZZ)1/3,  and  it  is  y  _  if  y  >  4(22)Va. 
e  cl  o  2  c2  o  2 

The  corresponding  values  of  r,  with  zq  vanishingly  small,  are 
obtained  by  subs titvtion  of  y^  and  y  ^  for  y  into  Eq.  (29).  There 
is  then  obtained: 


•'.♦‘(VY 


S  |  (22)1/2  , 


=  yo +  3(y, 


’  ■  4 )  v„  s  1  <22>1/a 


The  first  equation  above  is  identical  to  Eq.  (35),  which  was  derived 
on  the  basis  of  symmetrical  buckling.  Equation  (48)  is  the  result  of 
Fung  and  Kaplan  [5], 

Quasi-static  snap-through  for  this  problem  is  summarized 
in  Table  3.  It  may  be  inferred  from  these  results  that  the  appearance 
of  the  cross  sections  of  the  potential  surface  close  to  the  y-axis  is  that 
of  a  valley  for  y  <  4  (22)1/3  which  changes  to  a  ridge  near 

O  Ct 

y  =  y  7  if  y  >  ~  (22)1/3 . 

6  C  U  O  b 


5.  UPPER  BOUND  FOR  ELASTIC  ARCHES  WITH  INFINITESIMAL 
IMPERFECTION 


In  contrast  to  the  resulcs  for  the  perfect  elastic  arch,  for 

which  an  exact  closed-form  expression  (r  ,  Eq.  (32) )  for  the  critical 

cp 

load  can  be  determined  for  all  values  of  y  (y  >4  f^)1^3)  a*  which 

o  o  2 


snap-through  can  occur,  the  critical  load  for  an  infinitesimally 


imperfect  arch  cannot  be  represented  in  closed  form  for  y  >  5/(6)^3 

o 

It  has  been  shown  in  the  preceding  section  that  r  ,  Eq.  (43),  is  a 

cm 

lower  bound  and  generally  a  good  approximation  to  rc  for  y^  >  5/(6)^3 

An  upper  bound  for  r^  was  proposed  in  [2]  ,  and  a  summary  of  the 

results  obtained  there  is  as  follows. 

For  any  y  >  5/(6)1/3  consider  the  load  r  which  causes 
o  eg 

the  saddle  points  to  disappear ,  i.  e.  ,  the  load  which  causes  the  z  of 
Eq.  (39)  to  vanish.  It  is  found  that  this  load  is 


rcg  =  yo  +  3 'Vo8  -  4>1/a  •  <4’> 

Equation  (49)  has  the  same  right-hand  side  as  Eq.  (48),  but  the 

interpretation  here  is  different.  When  the  saddle  points  disappear, 

they  coalesce  with  one  of  the  critical  points  on  the  y-axis.  For 

y  <  —  (22)1/3 ,  this  critical  point  is  the  local  maximum,  but  for 
o  L 

y  >  —■  (22)1/3  f  the  critical  point  is  the  first  minimum  M.  .  When 
o  L  1 

this  occurs,  is  no  longer  a  stable  equilibrium  point  of  the  potential 

surface  with  y  and  z  coordinates,  and  dynamic  snap-through  must 
eventually  occur. 

A  critical  condition  will  occur  when  either  the  potential  U 


is  zero  at  the  local  maximum  or  becomes  unstable.  The  first 
condition  is  governed  by  r  and  the  second  by  rCg*  The  upper  bound 
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V  J  . 

I-  =&Zi 

4r-'- 


on  r  is  determined  by  whichever  of  these  two  values,  r  or  r  , 
c  7  cp  eg 

is  the  least  for  a  given  . 

These  results, and  the  other  results  presented  here  on 
dynamic  buckling  of  an  infinitesimally  imperfect  elastic  shallow  arch, 
can  be  summarized  as  follows: 


y  5  1.224 

o 


no  snap-through 


1.224  <  y  S  2.04 
o 


r  =  r 
c  cp 


2.  04  <  y  5  3.  19 
o 


r  5  r  5  r 
cm  c  cp 


3. 19  >  y 


r  —  r  5  r 
cm  c  eg 


A  comparison  between  this  upper  bound  and  the  calculated 
values  of  rc  is  shown  in  Table  1  and  graphically  in  Fig.  6.  The 
upper  bound  is  considerably  larger  than  r .  The  results  of  the  present 
investigation,  as  discussed  above  in  part  (b)  of  Section  4,  place  a  new 
interpretation  on  this  upper  bound,  however.  It  is  now  seen  that  the 
upper  bound  values  for  r^  in  (50)  may  be  those  actually  approached  by  some 
arches  with  very  small  viscosity  (x  =  0+)  while  the  lower  bound  values 
for  r^  in  (50)  are  those  actually  approached  for  elastic  or  zero  vis¬ 
cosity  (X  =  0)  arches.  In  both  cases  we  must  consider  that  the  arches 


are  slightly  imperfect  (z  =  0+).; 

o 


6.  ARCHES  WITH  FINITE  GEOMETRICAL  IMPERFECTION  (z  >  0+) 

o 

For  the  case  of  finite  geometrical  imperfections,  the 

critical  loads  can  be  found  from  the  numerical  solution  of  Eqs.  (17), 

(18)  and  (19)-:  The  critical  loads  for  z  =0.25  and  z  =  0  are  shown 

o  o 

as  dotted  curves  in  Fig.  5.  The  critical  loads  for  the  X  =  ®  curve  are 
obtained  from  the  maximum  point  of  the  r  -  y^  curve  evaluated  from 


Eq.  (22), 


As  compared  with  the  case  of  the  perfect  arch,  the  critical 


load  is  reduced  considerably  due  to  the  imperfections.  The  reduction 

increases  with  the  value  of  y  .  The  effect  of  finite  imperfections  for 

various  values  of  yQ  is  shown  in  Fig.  8.  As  expected,  the  critical 

load  decreases  as  the  value  of  initial  imperfection  increases.  In  accord 

with  our  previous  discussion,  a  discontinuity  in  the  critical  load  is 

found  at  z  =0  for  the  cases  y  =  3  and  4  but  not  for  y  =  2. 
o  7o  *o 

7.  PRINC IPA L  C PNC  LUSIONS 

In  summary,  we  have  observed  several  types  of  discontinuous 
behavior  of  the  calculated  snap-through  load  r  as  a  function  of  the 
imperfection  amplitude  z^  and  the  viscosity  parameter  x  .  This  be¬ 
havior  depends  also  upon  the  rise  y  ,  as  is  discussed  throughout  the 
text,  but  in  order  to  focus  on  the  essential  results,  we  will  discuss  the 
function  r^z^,  x)  within  the  range  of  yQ  upon  which  we  calculated. 


i 

r 

| 

•t 

f 

I 

a)  The  introduction  of  an  "infinitesimal"  imperfection  t 

*4*  ' 

causes  a  jump  or  discontinuous  drop  in  r  :  r  (0  ,  x)  <  r  (0,  x)  for  I 

C  C  C  r 

X  ^  0.  Larger  imperfections  lower  r^  continuously  (Fig.  8). 

b)  The  critical  loads  r  (0+,  0)  were  found  to  be  close  to 

c 

but  slightly  above  the  Hoff  and  Bruce  value  (Eq.  (43))  for  the  elastic 
arch  (see  Table  1  and  Fig.  6). 

c)  The  critical  loads  rc  (0+,  0+)  were  found  to  be  between 
the  upper  bound  loads  (Eq.  (50))  by  Simitses  and  the  lower  bound  loads 
(Eq.  (43))  by  Hoff  and  Bruce.  A  rapid  increase  of  the  critical  load  is 
observed  by  introducing  very  small  viscosity,  i.  e., 

rc(0+,  0+)  >  rc(0+,  0),  (see  Fig.  6). 

d)  For  finite  values  of  x  and  z  (Fig.  7),  r  increases  with 

o  °  c 

increasing  x  to  a  maximum  value  determined  by  the  quasi- static 
snap- through  load. 

e)  Computation  of  r^  for  both  zq  and  X  extremely  small 

can  lead  to  anomalous  results  for  r  because  the  arch  can 

c 

oscillate  for  a  long  time  before  snap- through. 

A  related  problem  of  dynamic  buckling  of  sinusoidal,  simply- 
supported  shallow  arches  has  been  considered  by  Lock  [8]  .  In  ht3 
analysis,  Lock  considers  what  is  called  in  the  present  paper  the 
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perfect  arch,  that  is  z  =0,  but  he  includes  infinitesimal  initial 

o 

antisymmetric  velocity  disturbances,  that  is  z^  of  the  order  of 
5  X  10  ,  which  are  not  considered  here.  These  antisymmetric 

initial  velocity  components  serve  in  Lock's  investigation  to  produce 
antisymmetric  snap-through.  Moreover,  in  Lock's  analysis,  energy 
dissipation  is  introduced  by  external  damping  rather  than  by  material 
viscosity  as  is  done  in  the  present  paper. 

Since  the  conditions  of  imperfection,  load  application  and 
damping  used  in  Lock's  work  and  the  present  work  are  different,  there 
is  no  reason  to  expect  the  same  conclusions.  The  following  points  can 
be  made  in  comparison  of  his  results  and  our  results: 

(1)  in  the  elastic  case,  Lock  finds  that  the  dynamic  buckling 

load  r  has  a  discontinuous  jump  at  y  =  4.375.  For  y  >  4.375, 
c  o  o 

the  value  of  r  can  be  higher  than  the  quasi-static  buckling  load 
c 

r  This  behavior  is  not  seen  in  the  present  case  for  y  s  4.5, 

eg  o 

(2)  In  the  case  of  small  damping.  Lock  finds  that  the  buckling 
load  for  y^  a  4.0  is  essentially  the  quasi-static  buckling  load.  The 
same  result  is  also  observed  in  our  calculation  (see  Fig.  6). 
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1J.  ABSTRACT  \ 

Dynamic  snap-through  or  dynamic  buckling  of  imperfect  viscoelastic 
shallow  arches  with  hinged  ends  is  considered  under  step  loads  of  infinite 
duration.  Attention  is  principally  devoted  to  the  influence  both  of  small  imperfect¬ 
ions  and  of  small  amounts  of  damping,  acting  together,  on  the  critical  loads.  For 
the  problem  considered,  the  Voigt  model  is  used  for  viscoelasticity,  the  deflection 
is  represented  by  the  first  two  harmonic  modes,  and  imperfections  have  the  shape 
of  the  second  (antisymmetric)  mode.  Results  obtained  by  numerical  integration  of 
the  differential  equations  show  that  the  critical  load  exhibits  a  jump  discontinuity 
in  the  limit  bot^h  for  vanishing  imperfection  and  for  vanishing  viscosity.  Critical 
loads  for  slight  imperfect  and  elastic  (inviscid)  arches  are  slightly  higher  than 
those  from  the  saddle  point  formula  of  Hoff  and  Bruce  (J.  Math.  Physics.  32.  1954, 
276),  confirming  that  the  formula  gives  a  lower  bound  on  the  critical  load. 
However,  critical  loads  for  arches  with  slight  imperfection  and  slight  viscosity 
are  considerably  higher  than  for  the  elastic  arches.  Another  closed-form  expres¬ 
sion  is  shown  to  be  in  good  agreement  with  these  results.  For  finite  amounts  of 
viscosity,  the  critical  toads  tend  rapidly  to  the  values  obtained  for  infinite 
viscosity,  which  are  the  same  as  the  critical  loads  for  quasi-static  buckling. 

Apart  from  the  jump  discontinuity  at  zero,  the  critical  load  for  any  viscosity 
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